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Abstract 

This paper compares two methods of on-orbit alignment of vector attitude 
sensors. The first method uses the angular difference between simultaneous 
measurements from two or more sensors. These angles are compared to the 
angular differences between the respective reference positions of the sensed 
objects. The alignments of the sensors are adjusted to minimize the 
difference between the two sets of angles. In the second method, the sensor 
alignment is part of a state vector that includes the attitude. The align- 
ments are adjusted along with the attitude to minimize all observation residu- 
als. It is shown that the latter method can result in much less alignment 
uncertainty when gyroscopes are used for attitude propagation during the 
alignment estimation. The additional information for this increased accuracy 
comes from knowledge of relative attitude obtained from the spacecraft gyro- 
scopes. This paper presents the theoretical calculations of this difference 
in accuracy. Also presented are numerical estimates of the alignment 
uncertainties of the fixed-head star trackers on the Extreme Ultraviolet 
Explorer spacecraft using both methods. 
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1. INTRODUCTION 


This paper compares two methods of determining the in-flight alignment estima- 
tion of vector-type attitude sensors. The two methods of alignment estimation 
will be referred to as attitude-independent and attitude-dependent and are 
outlined below. The estimated accuracy of the two approaches will be dis- 
cussed. 

Vector-type attitude sensors are those sensors whose output are vector mea- 
surements of the lines-of-sight to some reference object. Examples of this 
type of sensor are Sun sensors and star trackers. A minimum of two vector 
observations are needed for attitude determination. The analysis presented in 
this paper will be restricted to observations from two fixed-head star track- 
ers (FHSTs) for simplicity. 

Section 2 discusses the theoretical differences in the accuracy of these two 
approaches. Section 3 presents a numerical example of both methods using the 
Extreme Ultraviolet Explorer (EUVE) spacecraft as a typical mission case. 

ATTITUDE- INDEPENDENT ALIGNMENT ESTIMATION 

The attitude-independent method of sensor alignment estimation discussed here 
is based on an algorithm first presented in a paper by Shuster, Chitre, and 
Niebur (Reference 1) and later refined by Bierman and Shuster (Reference 2). 
This method of sensor alignment estimation uses the angle between two vector 
observations as its basic observation. This scalar observation is compared to 
the angle between the corresponding reference objects, and the sensor align- 
ment is adjusted to minimize the difference. The method requires simultaneous 
measurements in each sensor. No a priori knowledge of the attitude is used 
nor does the algorithm solve for the attitude. The algorithm seeks to mini- 
mize the overall deviation of the sensor alignments from their prelaunch 
values. 

ATTITUDE-DEPENDENT ALIGNMENT ESTIMATION 

In an attitude-dependent alignment estimation method, the spacecraft attitude 
is part of a state vector that includes the sensor alignments. The attitude 
is either solved-for along with the alignment or is treated as a known quan- 
tity. The primary observation quantity is the unit vector measurement from 
the sensors as opposed to the scalar measurement of the attitude-independent 
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method. Residuals are computed between estimated observations based on the 
state vector and the actual observations. The state vector is adjusted to 
minimize the residuals. Two procedures are in common use to obtain this 
minimization: the sequential filter and the batch filter. The sequential 
filter updates its estimation of the state vector at discrete times using a 
previously estimated state that has been propagated from the time of the last 
sensor observation. The batch filter considers a collection of sensor obser- 
vations in toto and seeks to minimize all residuals simultaneously. Both 
procedures require knowledge of the motion of the spacecraft between observa- 
tions. This knowledge is usually obtained from gyroscope measurements but 
could be inferred from a dynamics model of the spacecraft. The analysis in 
this paper will restrict itself to the batch filter and will assume dynamics 
information is available. 


2. THEORETICAL EVALUATION 

As stated in the introduction, this analysis will assume a spacecraft with two 
FHSTs. The state vector (References 3 through 5) to be estimated is 


x(t) 


A0(t) 

m 

l 


where A0(t) is a vector of attitude error angles around the spacecraft body 
— * 

axes and m^ (for i = 1 and 2) is a vector of sensor misalignment angles around 
the i th FHST axes (References 4 and 5). 

A batch least-squares estimate of the state at an epoch time t , ignoring the 
effects of dynamics noise and consider parameters, has the covariance 
(Reference 3) 

P(t ) s e[ x ( t ) x T (t ) 1 = W" 1 
o L 0 o J n 

where W is the normal matrix, 

n 


w = w + f t w f 

n o 
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This quantity contains attitude information, 
sensor alignments, B i and B 2> are given by 


For simplicity, assume the 


B 

2 


B i = * 3 x 3 * 3 x 3 identity matrix 


0 0-1 
0 1 0 
10 0 


- rotation by 90 degrees (deg) about the y-axis 


Thus, the sensor boresight, which is along the z— axis in the sensor frame, is 
along the spacecraft z- and x-axis for FHST-1 and FHST-2, respectively. The 
partial derivatives of the observations with respect to the state vector 
components are needed to compute the normal matrix. Again for simplicity, 
assume the observed stars are on the sensor boresights. The partial deriva- 
tive matrix for the scalar observation is 
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For the assumed geometry, the variance of the errors in the scalar observa- 
2 

tion, s, is 20-^. The partial derivative matrix for the vector observation is 
(References 4 and 5) 


ay 

M 

M 

o 1 

2x3 

ax 

MB 

2 

0 1 
2x3 

M 


where 


M = 


0 1 0 
-10 0 


The attitude and misalignment variances will be computed for four cases. For 
either a scalar or vector observation type, either a single measurement or 
else two measurements separated by a 90-deg attitude maneuver about the space- 
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The matrix W is the inverse of the a priori covariance of x and F T W F repre- 

O 

sents the information contained in the measurements. Assuming no a priori 
attitude knowledge and a priori knowledge of sensor misalignments with 
standard deviation <r per axis gives 


r o 

o 1 

3x3 

3x6 

0 

I 

L 3x6 

6x6 J 


where 0 denotes an n x m matrix of zeros and I is an n x n identity 

nxm nxn 

matrix. The form of the matrices F and W depends on the measurements. 

The i th FHST returns the two-component measurement 


Z 

i 


W 


w 


/ w 

/ w 


lz 

lz 


where W is the unit vector measurement of the star in FHST coordinates 

(References 4 and 5). The errors in the two components of Z^ are assumed to 

2 

be uncorrelated and to have equal measurement variances denoted by <r*. 

The observation at a given time processed by the estimator is a function of Z^ 

and Z . The attitude-independent observation is a scalar given by the inner 
2 

product of the two star vectors in the spacecraft body frame. 


where B and B are the sensor alignment matrices that rotate a vector from 
l 2 

the spacecraft body frame to FHST— 1 and FHST-2 frames, respectively. This 
quantity contains no attitude information. 

The attitude-dependent observation is a four-component vector 


y = 
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craft y-axis will be considered. The estimate epoch time will be taken to be 
the time of the first observation. The measurement weight matrix, W, will be 
assumed to be the inverse of the measurement covariance in all cases. The 
measurement error propagation matrix between the two observation times is 
given by (References 3 through 5) 


ft = 



0 


6x3 


| 0 1 

! 3x6 

I 

| I 

1 6x6 J 


where ft represents a 90-deg rotation about the y-axis. The four measurement 
0 

scenarios and their associated weight and total derivative matrices follow: 


1. A single scalar observation 



2. Two scalar observations separated by a 90-deg maneuver 


G 

s 

G ft 



I 

2x2 


3. A single four-vector observation 

F = G 


4. Two four-vector observations separated by a 90-deg maneuver 


F = 


G 

V 

G ft 

V 


W = 



I 

8x8 


For each scenario the normal matrix will be computed and then inverted to 
obtain the covariance matrix. In all cases a permutation of the rows and 
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columns of the normal matrix will bring it to block diagonal form, facili- 
tating matrix inversion. This permutation corresponds to rearranging the 
state vector as 


*■[ 


= A0 


lx 


A0 


m 


2x 


A0 


» » m , 

ly 2y 


1 z 


2z 


] T 


The normal matrix with the permuted rows and columns will similarly be dis- 
tinguished by a prime, W' . 

n 

SCENARIO 1: A SINGLE SCALAR OBSERVATION 

First consider scenario 1, a single scalar observation, for which 


f t w F = \ <r -2 g t g 

Z m 8 s 


o 

3x3 

0 

3x3 

0 1 

3x3 

1 -2 

= x <r 
Z m 

0 

3x3 | 

D 
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0 1 
3x3 ! 

-D 
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where D - diag(0, 1, 0) and diag(...) denotes a diagonal matrix with the 
indicated elements along the diagonal. The permuted normal matrix in this 
scenario is 


W' 

n 


-2 


<r 




o 


' 0 '_ 

j (a + 1/2) -1/2 

1 -1/2 (a + 1/2) i 


where 


a = (<r / <r ) 2 

m a 
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These variances, however, have the same limits for <r » <r and c « <r as in 

m a m a 

scenario 1. 

SCENARIO 3: A SINGLE FOUR-VECTOR OBSERVATION 
Scenario 3, a single four-vector observation, gives 


where 


- 2 

= < T 
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The permuted normal matrix is 


W' = cr 
n m 


-2 



5 a* 


l a ] 


Inverting this matrix gives the attitude variances 
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The variances of the state vector components are the diagonal elements of the 
inverse of W' . Thus, in this scenario the attitude error variances are all 

n 

infinite, which is to be expected with attitude- independent observations. The 

four misalignment variances, <r (m^), ^ ^ ^ m 2x^’ anc * ^ m 2z^’ are a ^ 

equal to their a priori values, <r 2 , which means that the single scalar obser- 
vation does not contain information about these misalignments. The y compo- 
nent of the misalignment variances for each of the two sensors is given by 

<r 2 + 1/2 cr 2 ) 

m a 

2 ~ 2 
cr + <r 


2 f . 2 

<r { m ) = <r 

iy a 


It is instructive to examine two limiting cases: the case of sensor data 

uncertainty much less than the a priori state vector uncertainty (o^ « <rj and 
the opposite case of poor measurement accuracy (cr » o^). The y-axis align- 
ment variance for scenario 1 has the limit for <r^ » and (1/2) <r^ for 
tr « <r The former result is reasonable since the measurements do not im- 

prove the a priori estimates in this limit. In the opposite limit of accurate 
measurements, the y-axis misalignments of the two FHSTs are in some sense 
averaged. These results are independent of the number of observations. 

SCENARIO 2 : TWO SCALAR OBSERVATIONS SEPARATED BY A MANEUVER 

For scenario 2, two scalar observations separated by a 90-deg maneuver, 


T 

F W F 


1 -2 ( „T 

= - <T G 

2 m ^ s 


G + $ T G T G * 

s s s 


) = <r - 2 


T 
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since it is easily seen that 


* T G T G * = g t g 

S3 S S 

Thus, the only change in W' from scenario 1 is to replace each 1/2 in the 
matrix by 1, and the only variances that are modified are 


2 f . 2 
o' (m ) = cr 
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<r 2 (A0 ) = <r 2 (A0 ) » <r 2 + <t 2 

X z a m 


<r (A0 ) 
y 





These results can be understood intuitively as follows. An observation of a 
star on the boresight of FHST-1, which is aligned with the spacecraft z-axis, 
provides information about the x- and y-axis attitude errors. Similarly, an 
observation of a star on the boresight of FHST— 2, which is aligned with the 
spacecraft x — axis, provides information about the y— and z— axis attitude 
errors. As there is twice as much information about the y — axis attitude, the 
variance of the y-axis attitude is half as large as that of the other two 
axes. The x- and z-axis uncertainties are the root-sum-square of the a priori 
alignment uncertainty and the measurement uncertainty. 

The alignment variances are the same as in scenario 1; therefore, no align- 
ment information is lost or gained by combining the four components of the 
vector observations into a single scalar observation for measurements taken at 
a single attitude. 

SCENARIO 4: TWO FOUR- VECTOR OBSERVATIONS SEPARATED BY A MANEUVER 

The final scenario, two four-vector observations separated by a 90-deg maneu- 
ver, has 


F T W F = <r -2 
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+ t G G $ 

V V 


where 


$ t g t g * = 

V V 


* T r $ 

0 00 0 

» T r 1 
0 e» 

T 

r $ 

L 0m 0 

r 

mm 


It is easy to see that 


$ T r * = r 

0 00 0 00 


however. 
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$ 


T 
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0m 


0 0 0 -1 0 

0 10 0 1 

-1 0 0 0 0 


0 


0 

0 


* r 


0m 


Thus, in contrast to the scalar observation case, the structure of the normal 
matrix will be different when a spacecraft maneuver is performed during align- 
ment. The normal matrix is 


W' 

n 


-2 


O 


m 


2 1 0-1 
1 (a + 2) -1 0 

0-1 2-1 
-1 0 -1 (a + 2) 


0 


j4 2 2 

\2 (a + 2) 0 

l 2 0 (a + 2) 


0 


lal 


5 CL i 


2 2 

It can be seen that <r (m^) = <r a in this scenario as in all the others. This 
result is simply that observations of stars on the boresight can never improve 
knowledge of the misalignment component about the boresight. The y-axis 
attitude and misalignment variances are the same as scenario 3 with <r 2 being 

2 m & 

replaced by (1/2) o reflecting the presence of twice as many measurements. 
Thus 


<r 2 (A0 ) 
y 





The y-axis misalignment uncertainties are the same as those in scenario 2, 

using scalar observations formed from the same measurements. The new feature 

of scenario 4 is the coupling of A0 , A0 , in , and m in the 4 x 4 block in 

x z lx 2x 

the upper-left corner of VT . Inverting this submatrix gives the variances 
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<r 2 (A0 ) = <r (A0 ) = O' 


+ 2 <r 


2 <r 2 + 2 <r 2 
m a 


o' (m ) = (<T + O' ) 

i x a ® 


The x-axis misalignment variance of the FHST is the harmonic mean of the 
measurement variance and the a priori alignment variance. It is, therefore, 
less than either variance and approaches the smaller of the two in the limit 
that the other becomes very large. The limits of the x- and z-axis attitude 
variances are (1/2) <r 2 for <r » <r and <r 2 for <r « <r . The latter is, like 

m nt a ro m a 

cr 2 (m ), finite as <r tends to infinity. Thus, unlike any other scenario 
i x & 

considered, the attitude and alignment knowledge obtained from the FHSTs is 
much better than the a priori knowledge of the misalignments. 

THEORETICAL EVALUATION SUMMARY 

Table 1 summarizes the alignment variance results for the four scenarios. As 
the table shows, the only significant difference in the alignment variance 
results between the attitude- independent and the attitude-dependent methods is 
due to the attitude maneuver. 


Table 1. Alignment Variance Results Summary 
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3. A NUMERICAL EXAMPLE 


This section presents a numerical example of the results obtained in Section 
2. The spacecraft that will be used in this example is EUVE, an astronomical 
satellite that is scheduled to launch in August 1991. Its primary attitude 
sensors are two FHSTs and three-axis gyroscopes. During the first 6 months, 
EUVE will perform a full-sky ultraviolet survey and will be spinning at rate 
of three revolutions per orbit around the spacecraft x-axis. After this 
survey phase, the spacecraft will be held at constant attitudes to measure the 
spectrum of various targets. The mission profile of EUVE, therefore, provides 
the opportunity of FHST alignment calibration under both maneuvering and 
stationary conditions. Although the alignment and attitude motion of EUVE 
differs from the analytical model presented in Section 2, they are sufficient- 
ly similar to substantiate the analytical results. 

NUMERICAL MODEL 

The two EUVE FHSTs are NASA standard star trackers with an 8-by-8 deg field- 
of-view (FOV). Their alignment is given as a 3-2-3 Euler rotation sequence. 
The rotation angles are given in Table 2. 

Table 2. Nominal FHST Alignments 


Sensor 

Rotation 

Angles 

(Degrees) 


0 

0 

0 


l 

2 

3 

FHST-l 

59.90 

105.60 

0.0 

FHST-2 

128.10 

105.60 

0.0 


The assumed a priori alignment uncertainties are 2.68 x 10~ 4 radians (55.3 
arc-seconds) (3<r) for each FHST axis. The measurement noise assumed is 
1.944 x 10 4 radians (40.1 arc-seconds) (3<r). In the attitude-dependent case, 
the a priori attitude uncertainties are assumed to be 1 deg (3<r) for each 
axis. The measurement rate for the EUVE FHSTs is one observation per 0.256 
seconds. To simplify the interpretation of the results, the gyroscopes are 
assumed to be perfect. 
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The attitude-independent alignment accuracy estimate is made using a program 
that implements the equations for the alignment covariance given in References 
1 and 2. The attitude-dependent alignment accuracy estimate is performed using 
the Attitude Determination Error Analysis System (ADEAS), as described in 
Reference 5. In both cases, a one-orbit simulation is run. Although the 
sensor coverage is not continuous because a realistic star catalog is used, 
the length of the simulation run and the data rate are sufficient to give 
several thousand observations. 

RESULTS 

The results for the four simulation runs corresponding to the four scenarios 
are given in Table 3. For simplicity, the arithmetic mean of the uncertain- 
ties for the two FHSTs are given. The actual differences are in all cases 
less than 1 percent. All values are three-sigma. 

Table 3. EUVE Alignment Uncertainties 


Axis 

Attitude- 
Independent , 
Nonrotating 
(Scenario 1) 

Attitude- 
Independent , 
Rotating 
(Scenario 2) 

Attitude- 
Dependent , 
Nonrotating 
(Scenario 3) 

Attitude- 
Dependent , 
Rotating 
(Scenario 4) 

X 

1.94 x 10" 4 
(40.0) 

1.90 x 10 -4 
(39.2) 

1.93 x 10' 4 
(39.8) 

0.32 x 10 -4 
(6.52) 

Y 

2.65 x 10" 4 
(54.7) 

2.03 x 10" 4 
(41.8) 

2.66 x 10" 4 
(54.8) 

1.50 X 10* 4 
(30.9) 

Z 

2.68 x 10" 4 
(55.3) 

2.07 x 10~ 4 
(42.6) 

2.68 x 10" 4 
(55.3) 

1.17 x 10" 4 
(24.1) 


NOTE: Values are in radians (arc-seconds). 


OBSERVATIONS 

Overall, the results in Table 3 and Section 2 agree. The smallest uncertain- 
ties are in the case of using an attitude— dependent alignment method during a 
maneuver. The differences between the results presented in Table 3 and the 
results from Section 2 are mainly due to the different geometry of the two 
cases. Figure 1 shows the fundamental geometry of the analytical and numeri- 
cal models. 
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Analytical model 


ROTATION 

AXIS 



Numerical model 


ROTATION 

AXIS 



FHST-1 FHST-2 


Figure 1. Geometry of the Analytical and Numerical Models 






In the analytical model, a rotation about the FHST y-axis is equivalent to a 
change in the angle between the sensors. The corresponding axis in the numer- 
ical model is the x-axis. Table 1 indicates that the uncertainty of this axis 
should be 1.90 x 10~ 4 radians (39.1 arc-seconds). This result is in very good 
agreement with the first three cases. Allowing for the difference in sensor 
coordinate definition, the uncertainties of the y- and z-axes for the two 
nonrotating cases are also in good agreement with Table 1 expected results. 
The smaller uncertainties in the rotating cases are due to the finite FHST FOV 
size in the numerical example. In the nonrotating cases, the FHST observes 
the same star for the entire simulation. Under this condition, the y and 
z-axes alignments are unobservable. This unobservability is also the conse- 
quence of assuming the star to be on the boresight as in the analytical eval- 
uation. In the rotating numerical examples, however, the realistic star 
catalog used results in star measurements over the whole FOV and the align- 
ments of the y- and z-axes are, therefore, observable. The variances of all 
three axes alignments in the attitude- independent case are larger than 
(1/2) cr 2 . 

a 

4. CONCLUSION 

The reason for the superior performance in the fourth measurement scenario of 
Section 2 is that the axis of the 90-deg attitude maneuver has been assumed to 
coincide exactly with the spacecraft y-axis. Since an actual estimator would 
obtain the angular rotation from a set of rate-integrating gyroscopes, this 
assumption is equivalent to the definition that the gyroscope axes are per- 
fectly aligned with the spacecraft body axes. An overall rotation of all the 
sensors, including unspecified payload instruments, is indistinguishable from 
a spacecraft attitude rotation. Because of this perfect correlation, if the 
sensors and payload instruments can be coaligned, the freedom exists to arbi- 
trarily choose one sensor to be perfectly aligned, which is equivalent to 
defining the spacecraft body axes in terms of this sensor’s axes. The 
attitude-dependent algorithm, as modeled in Section 2, implicitly defines the 
spacecraft axes in relation to the gyroscope axes, so the improved alignments 
of the FHSTs in scenario 4 are actually owing to their alignment with respect 
to the rotation vector as measured by the gyroscopes. 

An estimator that uses scalar observations as in Section 2 to align the atti- 
tude sensors does not include any attitude information. There is, therefore, 
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no measurement of the sensor alignment with respect to perfectly aligned 
gyroscopes. Because of the correlation between an overall sensor rotation and 
the attitude, there are three degrees of freedom in the sensor alignment that 
are unobservable. The resulting uncertainties in the alignment estimates, 
therefore, cannot be made arbitrarily small and will strongly depend on the 
prelaunch measurements. If the alignment uncertainty of one of the sensors is 
assumed to be zero, this sensor would become a reference sensor similar to the 
gyroscope in the attitude-dependent model. The alignment uncertainty of the 
other sensor should then be similar to the attitude-dependent results. To 
illustrate this effect, the scalar observation simulation software was exe- 
cuted assuming one of the sensors to be perfectly aligned. The results of 
this simulation are given in Table 4. Also given in Table 4 is the attitude- 
dependent, rotating case results from Section 3 for comparison. All values 
are three-sigma. 


Table 4. EUVE Alignment Uncertainties Showing 
Reference Sensor Case Results 


Axis 

Attitude- 
Dependent , 
Rotating 
(Scenario 4) 

Attitude- 
Independent, 
Rotating 
Reference Sensor 

X 

0.32 x 10" 4 

0.30 x 10“ 4 


(6.52) 

(6.19) 

Y 

1.50 x 10" 4 

1.32 x 10' 4 


(30.9) 

(27.2) 

Z 

1.17 x 10" 4 

1.50 x 10~ 4 


(24.1) 

(31.0) 


NOTE : Values are in radians (arc-seconds). 


It can be seen that the two cases are in good agreement. This result would 
indicate that the distinguishing factor for producing small alignment uncer- 
tainties is not attitude-dependence or -independence but whether there is a 
reference sensor that serves to define an on-orbit spacecraft coordinate 
frame. 

Whether, for a particular mission, an alignment algorithm that defines an 
on-orbit spacecraft frame or one that maintains the prelaunch frame should be 
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used depends on the type of payload. When the payload can align itself with 
respect to the attitude sensors, there is no need to maintain the prelaunch 
spacecraft reference. The ambiguity of the attitude reference can be removed 
by defining a new on-orbit reference, e. g. , the gyroscope axes, and much 
improved attitude accuracies can result. However, when the payload is atti- 
tude sensitive but not sufficiently so to allow coalignment, an algorithm that 
preserves the prelaunch frame is preferable. In this case, the alignment of 
the payload is known relative only to the prelaunch reference. A method that 
does not make use of a reference sensor minimizes the deviation of the on- 
orbit sensor alignment from the prelaunch measurement and, therefore, will 
result in the best estimate of the payload attitude. 
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